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Abstract 

The mechanisms of the magnetization switching of magnetic multilayers 
driven by a current are studied by including exchange interaction between 
local moments and spin accumulation of conduction electrons. It is found 
that this exchange interaction leads to two additional terms in the Landau- 
Lifshitz-Gilbert equation: an effective field and a spin torque. Both terms 
are proportional to the transverse spin accumulation and have comparable 
magnitudes. 

The concept of switching the orientation of a magnetic layer of a multilayered structure 
by the current perpendicular to the layers was introduced by Slonczewski [ffl] and Berger M , 
and has been followed up by Waintal et al. [§]. The central idea is that for a noncollinear 
configuration of the moments of the magnetic layer the current induces a torque acting on 
the spins of the conduction electrons which in turn transmit this torque to the background 
magnetization of the magnetic layers through the exchange interaction between conduction 



electrons and the local "d" electrons. An alternate mechanism of current induced switching 
was put forth by Heide et al @] in which the current across the magnetically inhomogeneous 
multilayer produces spin accumulation which establishes an energy preference for a parallel 
or antiparallel alignment of the moments of the magnetic layers; this magnetic "coupling" 
was posited to produce switching. Recent experiments have reliably demonstrated that 
the magnetization of a magnetic layered structure is indeed switched back and forth by 
an applied current f^Tj. However, it is unclear whether the magnetization switching is 
triggered by the current-driven effective field or by the spin torque mechanism or both. 

Here we examine the two views of current induced switching, spin torque and effective 
field, by solving the equations of motion for the spin accumulation and the local magne- 
tization. We find the two mechanisms do coexist; albeit in form very different from that 
envisaged by the above referenced authors. The salient difference between our treatment of 
spin diffusion and previous treatments P-jTH, lies in the inclusion of the exchange interac- 
tion between the spin accumulation and the magnetic background. With our results, we can 
understand these two mechanisms on an equal footing: both are simultaneously derived and 
both depend on the same set of parameters used for understanding the giant magnetoresis- 
tance when the current is perpendicular to the plane of the layers (CPP). Furthermore, we 
have introduced a new length scale for the transverse spin accumulation and clarified the 
ferromagnetic layer thickness dependence of the switching dynamics. 

Let us consider a magnetic multilayer with the current perpendicular to the plane of the 
layer (defined as x-direction) . The linear response of the current to the electrical field can 
be written as a spinor form, 

j(x) = CE(x) - D^-, (1) 

where E(x) is the electric field, j, C, D, and n are the 2x2 matrices representing the current, 
the conductivity, the diffusion constant, and the accumulation at a given position. The 
diffusion constant and the conductivity are related via the Einstein relation C = e 2 N(ep)D 
for a degenerate metal, where N(ep) is the density of states at the Fermi level. In general, 
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one can express these matrices in terms of the Pauli spin matrix cr, 

c = c i + <t ■ c, 



(2) 



D = D I + a ■ D, (3) 

and 

n = n I + <r ■ m, (4) 

where 2no is the charge accumulation and m is the spin accumulation. By placing Eqs. (2)- 
(4) into (1), we rewrite the linear response in terms of the electric current j e and magneti- 
zation current j m as 

j e = Re(Trj) = 2C E(x) - 2D ^ - 2D ~, (5) 

and 

j m = ReTr(crj) = 2CE(x) - 2D^ - 2D ^. (6) 

It is noted that we have chosen the unit e = fiB — 1 for the notation convenience so that 
the electrical current and the magnetization current have the same unit. For a transition 
ferromagnet, one defines the spin polarization parameter (3 as C = f3CoNld, where is the 
unit vector to represent the direction of the local magnetization. Similarly, we can introduce 
a spin polarization (3' for the diffusion constant D = /3'DoM.d- These two polarization 
parameters are not necessarily the same, i.e., when the density of states are different for 
spin up and down electrons, j3 ^ f3' . By inserting these relations into Eqs. (5) and (6), and 
eliminating the electric field and charge density, we obtain 



(3 Je M d - 2D 



(7) 



ox ox 

where we have dropped an uninteresting term proportional to the derivative of the charge 



accumulation dn /dx. This result is similar to that obtained by Heide [11] 



We now describe the equations of motion for the spin accumulation and local magnetiza- 
tion when we turn on the interaction between the spin accumulation and the local moment 
via the sd contact interaction, 

H mt = - Jm ■ M d . (8) 

With this interaction, the equation of motion for the spin accumulation is 

dm / t / -j- \ ^ r m 

— + J/fi m x M d = , (9) 

dt r s f 

where r s f is the spin-flip relaxation time of the conduction electron. The second term 
on the left hand side represents the processional motion of the accumulation due to the 
sd interaction when the magnetization directions of the spin accumulation and the local 
moments are not parallel. Since the conduction electrons carry a spin current given by 
Eq. (7), we replace by ^ + By using Eq. (7), we find 

1 dm d 2 m /, , <9 2 m\ m m x M_ d . . 

" — - M ^ "Sir " T2 \2 ' 



2D dt dx 2 ^ a V dx 2 J X 2 f Xj 



where we have defined A s / = J2D T s f and Xj = J2hD / J (T^]. The latter gives rise to a 



new length scale which governs the spin torque created by the current. The significance of 
this new length scale will be discussed later. 

The equation of motion for the local magnetization is 

— — = -7 M d x (H e + Jm) + aM d x — — , 11 
dt dt 

where 70 is the gyromagnetic ratio, H e is the magnetic field including the contributions from 
the external field, anisotropy and magnetostatic field, the additional effective field Jm is 
due to coupling between the local moment and the spin accumulation, and the last term is 
the Gilbert damping term. 

To solve for the dynamics of the spin accumulation and the local moment, we need to 
simultaneously determine them using Eqs. (10) and (11). The time scales are very different 
for the spin accumulation and the local moments. The characteristic time scales of the 
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former are of the order of r s / and h/J, see Eq. (9), i.e., of the order of picoseconds (10" 
seconds). For the local moment, the time scale is 7 c j~ 1 (if e + Jm±)~ 1 . For a magnetic field of 
the order of 0.1 Tesla, this time scale is of the order of nanoseconds. Therefore, as long as 
one is interested in the magnetization process of the local moments, one can always treat the 
spin accumulation in the limit of long times. The two dynamic equations are then simply 
decoupled: we first solves Eq. (10) with fixed local moments (independent of time) and set 
the left hand side of Eq. (10) to zero. Once the spin accumulation is obtained, we substitute 
it into Eq. (11) to solve the dynamics of the local moments. 

Before we proceed to solve for the stationary solution of Eq. (10), let us first discuss the 
general features derived from Eq. (10). We separate the spin accumulation into longitudinal 
(parallel to the local moment) and transverse (perpendicular to the local moment) modes. 
Equation (10) can now be written as 



cr my m| 

dx 2 A 2 



0, (12) 



sdl 



where X sd i = ^1 - ]3f3'\ s f, and 

<9 2 m_i_ m_i_ m_i_ x M a 



«9* 2 \l f Xj 



0. (13) 



The longitudinal spin accumulation my decays at the length scale of the spin diffusion length 
X s di while the transverse spin accumulation m_|_ decays as Aj. The spin diffusion length A s <« 
has been measured to be about 60nm in Co fll3| . We estimate Aj by taking the typical 
diffusion constant of a metal to be 10 -3 (m 2 /s) and J = 0.1 — 0.4 (eV) so that Aj is about 
1.2 nm to 2.4 nm. Thus, the transverse spin accumulation has a much shorter length scale 
compared to the longitudinal one. 

Before we apply Eqs. (12) and (13) to a multilayer structure, we take a look at the influ- 
ence of the spin accumulation on the local moment. As seen from Eq. (11), the longitudinal 
spin accumulation has no effect on the local moment. We may re- write Eq. (11) in terms of 
the transverse spin accumulation only, 

^ = - 7o M d x (H e + Jmi) + aM d x (14) 



To discuss the transverse accumulation we introduce a vector A such that Jm± = A x M^. 
If one considers a system with two ferromagnetic layers whose magnetization directions are 
not parallel to each other, the spin accumulation at one layer depends on the orientation 
of the other. Let us suppose that the above equation is used for the layer Fl, i.e., denote 
M d = M^. The magnetization of the other layer is labeled as M^. Without loss of 
generality, we can write the two components of the accumulation in the plane transverse to 
mJ } as Jm_L = aM^ ] x + b(M { j ] x mJ } ) x M^, where a and b are determined by 
geometric details of the multilayer. Placing this into Eq. (14), we find 

( ' M[! ' -loM^ x (H e + 6M< 2 >) - 7o aM« x ( M f x M«) + aM? x (15) 



dt ' U a " a j a v a a , - a ^ 

Thus the transverse spin accumulation produces two effects simultaneously ( one can call 

(2) 

them either fields or torques): one is bM^' the "effective field" which gives rise to a pre- 
cessional motion and the other is aM^ 1 ' x (Mjj x M^) which is called the "spin torque". 
Both terms lead to significant corrections to the original Landau-Lifshitz-Gilbert equation. 



It has been shown that both terms are capable to switch the magnetic moments ||14|| . Note 
the effective field introduced here looks as if it arises from the current induced coupling 
named NEXI, however it is different as NEXI was attributed to the longitudinal component 



of the spin accumulation ||11|| . In contrast we have shown that only the transverse spin 
accumulation must be taken into account and that the longitudinal accumulation does not 
produce any effect on the motion of local moments. An even more striking difference is 
Heide's finding that "the presence of a second ferromagnetic layer is not necessary" . This 
is because his longitudinal accumulation exists for a single F layer, while a second F layer 
with tilted magnetization is required for transverse accumulation and for our mechanism. It 
is notable that the "torque" term, first introduced by J. Slonczewski, appears on an equal 

(2) 

footing with the effective field as both are related to the transverse spin accumulation. 

We now explicitly verify that the solution of the transverse accumulation mi indeed has 
our proposed general form and we quantitatively determine the magnitude of the effective 
field (proportional to b term) and the " spin torque" (proportional to a term) entering 
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Eq. (15). To obtain a physically transparent solution of the spin accumulation, we choose 
an oversimplified case to perform our calculation so that the effective field and spin torque 
can be analytically derived. We consider a system consisting of a very thick ferromagnetic 
layer which is assumed to be pinned, a spacer layer which is infinitely thin so that the spin 
current is conserved across the layer when there is no spin flip scattering in this region, 
and a thin ferromagnetic layer backed by an ideal paramagnetic layer. In addition we make 
our calculation simpler by neglecting spin- dependent reflection at the interfaces. In such 
a system, we look for the solution of the spin accumulation in the thin Fl layer whose 
magnetization direction is at the positive z-direction. The magnetization direction of the 
pinned layer is = cos^e^ — sin^e^, where 9 is the angle between and = e z . 
From Eqs. (12) and (13), and by assuming the same X s di for the thin magnetic layer, Fl, as 
for the non-magnetic layer which backs it, we write the solution for the Fl layer as 

m z (x) =G 1 exp(-x/X sdl ) (16) 
m x(x) = G 2 exp (— x/l + ) + G 3 exp (—#//_) (17) 
m y (x) = —iG<i exp (—x/l+) + iG% exp (— x/lJ) (18) 



where L- 1 = , -rh- ± To determine the constants of integration, we assume the thick 
magnetic layer F2 is half metallic so that the current is fully spin polarized and we demand 
that the spin current is continuous across F2/N/F1 interface Jl5|]; we find 



/3j e - 2A>(1 - (3(3') ( ) = j e cos9, (19) 



A. 



di 



J+ l- 

and 



2^o + = 0, (20) 



2D (-z) + j*)= -Jesin9. (21) 



Thus we determine the constants to be 



Go 



and 



jeKcojP ~ COS 9 

' 2D (l-P/3') 
j e l + sin 9 

= AtDn ' 



j e /_ sin 6 1 
AiD n 



(22) 



(23) 



(24) 



Therefore, we find the transverse spin accumulation 



mi = — 



( Je 



\2D 



lm{l + e- x,l +)mf + Re(/+e- x// +)Mi 2) x M. 



(2) 



r(l)^ 



xM« 



(25) 



where we have used -sin8e x = mJ } x mJ } and sin 6^ = (M^ 2) x M^) x mJ J . We 
immediately see that the form of the spin accumulation given above is precisely the form we 
used in deriving Eq. (15). To obtain the coefficients a and b entering Eq. (15) we average 
this spin accumulation over < x < tp where tp is the thickness of the Fl layer and find 

Jje 



2D t F 



Im[^(l - e -*W+)] 



(26) 



and 



Jje 



■Re[Zi(l -e 



(27) 



2D t F 

It is noted that both a and b change sign under time reversal. The former agrees with 
that found in |1|-[3|]; while the latter has not been considered by these authors. To estimate 
a and b, we take the limit A s / ^> \j in which case 1+ — (1 + i)Xj/y2. By placing this into 
Eqs. (26) and (27), we find 

foje a Q ( 1 — C °S £ e ~^ 



(28) 



and 



hj e al ( sin£e 



y/2efj, B K 



(29) 



where £ = tp/(\/2\j), ao is the lattice constant, and we have reinserted the electric charge 
and Bohr magneton so that a and b have units of a magnetic field. If we take Xj = 20 A, 
a Q = 2A, j e = 10 11 A/m 2 , we find a = —1056 (Oe) and b = 457 (Oe) for a typical experiment 
with t F = 25A. 

In conclusion we have found that by considering the exchange forces between the con- 
duction electron spin and the background magnetization for the spin current perpendicular 
to the layers of a magnetic multilayer there exists the effective field and torque, both of which 
contribute to current driven reversal of the magnetization. We treat both terms on an equal 
footing and demonstrate that they have a common origin. Our solution differs in two impor- 
tant aspects from previous work: we find the longitudinal spin accumulation does not play a 
role in the switching, and the spin torque, as well as the effective field, arises from a region 
in the magnetic layer within ~ Aj of the interface. Therefore, the decay length in our theory 
is related neither to the phase of the wavefunction , nor to the spin diffusion length as 
in the effective field concept of switching || . We would like to acknowledge our very fruitful 
conversations with Yaroslaw Bazaliy, Piet Brouwer, Carsten Heide, Henri Jaffres, Barbara 
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